Mechanisms of 


. N82 23425 

m 


droplet combustion 
C. K. Law 

Department of Mechanical and Nuclear Engineering, Northwestern University 

Evanston, Illinois 60201 


Abstract 


The fundamental physico-chemical mechanisms governing droplet vaporization and combus- 
tion are discussed. Specific topics include governing .'quations and simplifications, the 
classical d^-Law solution and its subsequent modification, finite-rate kinetics and the 
flame structure, droplet dynamics, near- and super-critical combustion, combustion of mul- 
ticomponent fuel blends/emulsions/suspensions, and droplet interaction. Potential research 
topics are suggested. 


1. Introduction 


Research on the gasification, oxidation, and dynamics of fuel droplets commands both 
practical and fundamental interest to energy and combustion science. On the practical 
aspect one recognizes that petrole'>m oil constitutes a significant share of the world 
energy supply. Since these fuel oils are usually introduced into the combustor as sprays 
of droplets, it is reasonable to expect that the collective gasification c ' individual 
droplets would intimately influence the bulk spray vaporization and oxidat. n characteris- 
tics, which in turn determines the combustor performance. On the fundamental aspect drop- 
let combustion is a problem involving complex chemically-reacting multicomponent two-phase 
flows with phase change, rich in physical and chemical phenomena tynically of interest to 
the study of aerothermochemistry . 

Figures la to Ic show some possible droplet combustion modes. In Fig. la the fuel drop- 
let is motionless in a stagnant, gravity-free, oxidizing environment of infinite extent. 

The lack of either forced or natural convection implies that spherical symmetry exists. 

The basic mechanisms leading to the complete gasification of the droplet are heat and mass 
diffusion due to the existence of temperature and concentration gradients, the radial con- 
vection because of the continuous transfer of mass from the droplet surface to the ambi- 
ance, and chemical reaction in the flame region. The burning is of the diffusion-flame 
type in which the outwardly-diffusing fuel vapor and the inwardly-diffusing oxidizer gas 
approach the reaction zone in approximate stoichiometric proportion. Reaction between them 
is rapid and intense, implying that the reaction zone is thin and very little reactants can 
lea . through the flame. The chemical heat generated at the flame is transported both out- 
ward and inward to heat up the approaching oxidizer and fuel gases in order to achieve 
ignition. The rest of the inwardly-transported heat is used for droplet heating and to 
effect liquid gasification. 

In the presence of either forced and/or natural convection, a non-radial relative veloc- 
ity exists between the droplet and the surrounding gas. The shear stress exerted by the 
gas flow on the surface induces a recirculatory motion within the droplet (Fig. lb). For 
higher rates of external convection, flow separation occurs close to the rear stagnation 
point, creating wake regions both inside and outside of the droplet (Fig. Ic). The pres- 
_nce of the non-radial convection generally enhances the transport rates and thereby the 
gasification rate, although for sufficiently intense flows the envelop diffusion flame can 
be extinguished, leading to significant reductions in the gasification rate. It is clear 
that fluid motion can have strong influence on the combustion process. 


The bulk parameters of interest in droplet combustion studies include the droplet gasi- 
fication rate, the flame location and temperature, the droplet drag, ignition and extinc- 
tion limits, and the extent of pollutant formation such as that of NO^. The droplet 
diameters of interest are typically between 10 um and 100 ym; gasification of smaller drop- 
lets can be considered to be instantaneous during the period between spray injection and 
active combustion, while large droplets either tend to breakup or cannot achieve complete 
gasification during the available residence time. Droplet breakup also limits its maximum 
Reynolds number, Re, to about 100 for typical surface tension values of hydrocarbon fuels. 
Re being defined as the ratio of the characteristic inertia to viscous fo»"ces. It may also 
be noted that it is frequently desirable for the droplets within a spray to have a distri- 
bution of sizes and velocities in order to achieve an optimum spatial distribution through 
penetration as well as a controlled rate of gasification and thereby chemical heat release. 


The basic spherically-symmetric, pure-component, droplet combustion model was formulated 
in the fifties. This model has since been termed the d^-Law because it predicts that the 
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square of the droplet diameter decreases linearly with time. In subsequent years the basic 
understanding of this model has been extensively revised while other aspects of droplet 
combustion have also been investigated. Extensive reviews of these developments can be 
found in Refs. 1 to 3. 

The purposes of the present paper are to familiarize the reader with the influence and 
modeling techniques of chemical reactions in droplet processes, to briefly review current 
understanding of the basic physico-chemical mechanisms governing droplet combustion, and to 
suggest potential areas of research. In the next section some frequently-invoked assump- 
tions during droplet combustion modeling are stated and a set of simplified governing equa- 
tions presented. The d^-Law is then derived, and its limitations and subsequent modifica- 
tions in light of experimental results discussed. This is followed by special topics on 
finite-rate kinetics, droplet dynamics, high-pressure combustion, multicomponent droplet 
combustion, and droplet interaction. In Section 9 potential areas of research are 
suggested . 


2 . Assumptions and simplified governing equations 

The general equations governing droplet combustion are those of heat, mass, and momentum 
conservation, in both the gas and liquid phases, and their coupling at the interface. 
Auxiliary relations needed are the equations of state and those describing chemical reac- 
tions and interfacial phase change. The chemical reactions of interest are usually those 
of fuel vapor decomposition and oxidation in the gas phase. For very high boiling point 
fuels, their pyrolysis in the liquid phase may also be of relevance. These equations, for 
general combustion problems, can be found in the text by Williams.** 

We shall instead study in the following a simplified system which illustrates the basic 
diffusive-convective-reactive nature of droplet combustion. The major assumptions are as 
follows. 

(Al) Diffusion being rate-limiting, although second-order diffusion, that is the Soret and 
Dufour effects, are neglected. 

(A2) Isobaric processes because of the low Mach number flow. 

(A3) Gas-phase quasi-steadiness. Because of the significant density disparity between 

liquid and gas, the liquid possesses much inertia such that its properties at the sur- 
face, for example the regression rate, temperature, and species concentrations change 
at rates much slower than those of the gas-phase transport processes. Thus the gas- 
phase processes can be treated as steady, with the surface boundary variations occur- 
ing at longer time scales. 

(A4) Constant gas-phase transport properties. A single binary diffusion coefficient, 6, is 
used for all pairs of species. The specific heat Cp, the thermal conductivity coef- 
ficient X, and the product p£, are assumed to be constants, where p is the density. 

With the above assumptions, the gas-phase governing equations arc 


Continuity 

V^PU) ■ 0 

(1) 

Species 

v*(puYi - pfiVYi) « Wi , i>l,2,...N 

(2) 

Energy 

V*(puCpT - X7T) ■ q 

(3) 

State 

p » p(RO/W)T 

(4) 


where T is the temperature, p the pressure, w an average molecular weight, u the velocity, 
Yi the mass fraction of species i, w^ the chemical production rate of i, q the chemical 
heat generation rate, and N the total number of species. 

For each of the equations represented by Eqs. (2) and (3), the first and second terms on 
the LHS respectively represent convective and diffusive transport, while the RHS is the 
chemical source/sink term. From chemical kinetics it is well known that the oxidation of a 
fuel species leading to the formation of products occurs through a complicated sequence of 
intermediate reactions with a myriad of intermediate species. Thus, for example, the term 
Wj represents the sum of the mass rates of production and destruction, per unit volume, 
of species i from all of the reaction steps. 

The adoption of such a detailed kinetic scheme yields immensely complicated expressions 
for wj and invariably necessitates the use of numerical solutions. Furthermore, except 
for the simplest hydrocarbon fuels, the detailed oxidation kinetics of the heavier ones are 
largely unknown. Thus frequently it is both necessary and useful to approximate the 
detailed kinetics as a one-step overall reaction between the fuel, F, and oxidizer, 0, as 
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(5) 


Oi 




k 


0 

(PI + ootOl * Products 


where on is the stoichiometric mass coefficient of the oxidizer per unit mass of fuel 
consumed, and Q is the associated heat release. The consumption rate of P can then be 
expressed as 


■ - Ap 


0 F 


exp ( 


R°T 


( 6 ) 


which shows that the overall reaction rate is proportional to the concentrations of the 
reactants and is dependent on the temperature through the Arrhenius factor characterized by 
an activation energy E. The exponents a and b are the empirically determined reaction 
orders with respect to the fuel and oxidizer, while A is the frequency factor which is a 
measure of the collision frequency between the reactants. Since pi^xY^p, Eq. (6) can 
also be expressed as 

w - - Y* Y^’ exp {- — ) (7) 

P 0 F pO^ 

Finally, recognizing that in any given reaction the rate of chemical heat release as 
well as the rates of creation and destruction of all species must be stoichiometrically 
related, we have 


°0 '*F 


and 


( 8 ) 


<1 - - OWp 

for the reaction given by Eq. (5). 

Equations (2) and (3) are now well defined. It is clear that the highly non-linear 
nature of the reaction rate term, which also strongly couples the Y^ and T variations, 
introduces challenging mathematical difficulties and interesting physical phenomena to the 
problem. 

Considerable simplification in the analysis results if one adopts the Shvab-Zeldovich 
Formulation, ** which requires 

(A5) Unity Lewis number, or Le » X/Cpp« - 1. 

Thus a linear combination of Eqs. (2) and (3) eliminates the reaction term from all but one 
of the (N-t-l) equations. In particular, if we define a coupling function 



(9) 


where op ■ 1, then Eq. (2) can be replaced by 

Tj^-Conservation 7*(pur^ - p«vr^) » 0 i»0,F (10) 

which shows that is a conserved scalar during the chemical reaction of Eq. (5). 
Physically represents the sum of the thermal and chemical energies which remains 
constant when heat and mass are transported at equal rates. 

The problem is thus simplified to the solution of the cheraically-inert Eq. (10) together 
with Eq. (3). Sometimes it is possible to first solve for such that Y^, and thereby 
Eq. (3), become only a function of T. The problem is decoupled. 

To demonstrate this methodology and be more specific in our discussions, let us further 
assume 

(A6) Spherical symmetry, implying the absence of natural or forced convection. 

Then the governing equations become 


f**' 
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(r^ ai) 

- Y« Y^’ exp 1- E- 


- 

(13) 
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Vhere r is the radial co-ordinate, s indicates the surface, m»( pur^ )/( p6rs) , r^r/rgt 
T-CpT/0, 


D . (A (I£!I) 

p« qrO 


CpP'* \ 1 I E 


} exp (- 


-) 


R°Tr 


is the ratio of the characteristic flow time to the characteristic reaction time, and 
is a characteristic temperature in the reaction zone. 

Equation (11) shows that 

m » constant 

implying that the mass flow rate is a constant as it should be. 


(14) 


Equations 

(12) and (13) 

are to be solved subject to the boundary conditions 


At r = » ! 

1 

o 

>4 

H 

o 

>4 


(15) 


Yp - 0 


(16) 


T = T» 


(17) 

At r » 1 : 

dYo/dr » 

mYO 

(18) 


dYp/dr = 

-m(l-Yp) 

-19) 


dT/dr = 

litL/O 

(20) 


*1* Z 

Ts 

(21) 


where L is the latent heat of gasification. Equations (18) and (20) respectively state 
that oxidizer does not penetrate into the surface, and that all the heat conducted to the 
droplet is used for gasification. Equations (15) to (21) constitute seven boundary 
conditions to solve for the three second order differential equations, Eqs. (12) and (13), 
together with the burning rate parameter m. 

Solving Eq. (12), we have 

(1° + t) • [l2!!)e""’^'^ - 8(l-e“'"^'^) + T (22) 

OQ OQ • 

and 

(Yp+T) - (1-8) (l-e-'"/f) + T« (23) 

where 


8 


(T.-Tg) + L 
__ 


Thus the problem is reduced to the solution of Eq. (13) for T, with Yq and Yp gi.».. by 
Eqs. (22) and (23) . 




’V»v 


o 


3 . The d^-Law and its modifications 




A particularly simple solution exists if we make the following assumption. 

(A7) Reaction proceeds with an infinitely fast rate at an inf initesimally-thin flame 
located at ff. Implying that both reactants are totally consumed there. 

Thus we have 


Yq ■ Yp ■ 0, T ■ Tf , at r ■ rf 


( 24 ) 


Furthermore, since there is now no leakage of the reactants through the flame, the oxidizer 
concentration is identically zero in the inner region to the flame. Thus, in particular, 
at the droplet surface we have 

Yq ■ 0 at r » 1 (25) 


Applying Eqs. (24) and (25) to Eqs. (22) and (23), we obtain three equations from which the 
three parameters of practical interest, namely the burning rate m, the flamefront standoff 
ratio rf, and the flame temperature tf, can be solved as 

m « tn(l+B) (26) 


f tn ( 1 +Yq^ 0 q) 


and 


Cr,(Tf-T„) 

(-E — i } + {L + C (T -T ) } = Q 


where B « 


Cp(T.-Tg) + (V/Oo’Q 


(27) 


(28) 


(29) 


is a transfer number representing the ratio of the "driving force" for vaporization to the 
"resistance" to vaporization. The chemical source term, (Yo»/oo)Q» is usually much 
greater than the thermal term, Cp(T»-Tg). Equations (26) to (28) show that m, ff, 
and Tf are constants for a given fuel-oxidizer system. 

Finally, further assume that 

(A8) Fuel vaporization rate at surface is equal to fuel consumption rate at flame m»4itpur^. 
Then 


- 1_ (1 , p d^) 

dt 3 t s 

which is simply 

d(d^) 

1- - - K 

dt 


with 


(30) 


(31) 


K - tn(UB) (32) 

where p^ is the liquid density. The parameter K is known as the evaporation or burning 
rate constant. 


Integrating Eq. (31) yields 


so 


- Kt 


which forms the basis of the d -Law. 


( 33 ) 
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0, we obtain the total burning time as 


Setting ds ■ 



K 


Equation (34) shows that t varies quadratically with the initial droplet size, hence demon- 
strating the basic principle of atomization in order to achieve rapid gasification. 

It is also of interest to note that in Eq. (28) we have expressed Tf in such a way as 
to illustrate the fact that it is just the adiabatic flame temperature for the stoichio- 
metric reaction between the liquid fuel and the oxidizer-inert mixture. That is, the 
amount of heat release, Q, by burning unit mass of fuel is used (a) to gasify the fuel 
and heat the fuel vapor from Tg to Tf, and (b) to heat (Yo./oo)~^ amount of the 
oxidizer mixture from T« to Tf. 

Figure 2a shows the d^-Law temperature and concentration profiles, which are typical of 
diffusional burning. 

Because of its simplicity, the d^-Law has been extensively used in modeling droplet and 
spray combustion. In particular, the qualitative behavior of Eq. (33) has been found to be 
largely correct while the burning rate constant can generally be predicted to within a 
factor of two, allowing for the uncertainty in selecting transport property values. 

However, experiments®*^ on spherically-symmetric droplet combustion reveal certain 
qualitatively different behaviors, indicating the existence of serious wea)cnesses in the 
d^-Law. Figure 3 shows that, immediately after ignition, there exists a short period 
during which dg^ practically does not change. Furthermore, rather than a constant, ff 
increases significantly after ignition. In a Iow-Yq- environment the increase persists 
until burnout, while in a high-YQa environment rf levels off toward the latter part of 
the droplet lifetime. 


Both of these deviations are consequences of initial conditions which are not included 
in the d^-Law. The initial slow rate of surface regression is primarily due to the need to 
heat up the initially cold droplet to a temperature hot enough to sustain steady burning; 
frequently this temperature is close to the liquid's boiling point under the prevailing 
pressure. Du'ing this period most of the heat generated is used for droplet heating 
instead of gasification, thereby resultinj in a slow surface regression rate. The droplet 
also expands as it is heated up, hence compensating for the reduction in the droplet size 
from gasification. 


Since active droplet heating and intense gasification are somewhat mutually exclusive, 
droplet heating is mostly over in about the first 10% of the droplet lifetime. For the 
remaining period dg should vary approximately linearly with time. During multicomponent 
droplet combustion or high-pressure combustion the droplet heating period can be considera- 
bly longer, as will be discussed later. 


To model droplet heating one may simply include in the latent heat of vaporization, L, 
an additional amount representing the heat conducted into the droplet interior. This 
amount can then be determined by solving the temperature distribution within the drop- 
let.^*® 


)■ 




The observed behavior of ff can be explained by the fact that the fuel 'ipor present 
in the region between the droplet and the flame continuously change as the droplet and 
flame sizes change. The significance of this variation can be appreciated by considering 
a droplet immediateljr after it is ignited. At that instant the amount of fuel vapor pres- 
ent in the inner region should be of the same order as the amount present in the droplet 
vicinity before ignition. Since this amount is very small because of the low droplet 
temperature, the flame initially must lie close to the droplet surface. With subsequent 
gasification the fuel vapor present in the inner region increases and larger flames can be 
supported. Therefore only part of the fuel vaporized during this period is consumed at the 
flame, the rest is being accumulated in the inner region as the flame expands. This amount 
may also become depleted toward the latter part of the droplet lifetime. 

It is reasonable to expect that the accumulated amount is significant because although 
the gas density is low compared with the liquid density, the flame size can be substantial 
such that volume effect dominates. Furthermore the fuel vapor is accumulated at the 
expense of the finite droplet mass. Indeed it is easy to demonstrate® that the amount of 
the fuel vapor present, as given by results of the d ^-Law, is of the same order as the 
droplet mass. 

The existence of this accumulation process implies that mass conservation is actually 
violated in d^-Law. Rather, overall mass conservation for the fuel vapor should read. 
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Gasification rate at droplet surface 

- Consumption rate at flame 

+ Accumulation/Depletion rate in the inner region (35) 

The last term is absent in the d^-Law. 

When the accumulation process is included in the formulation^, then the experimentally 
observed behavior of r^ is predicted. The result that ff does approach a constant 
value for burning in high Yq« environments can also be explained by its relatively 
smaller flame size, which requries less fuel vapo.' for accumulation. 

An Important practical implication of fuel vapor accumulation is that since the fuel 
gasification rate is not equal to the fuel consumption rate, adoption of the d^-Law in 
spray modeling may result in grossly erroneous estimates of the bulk chemical heat release 
rate. In particular, in a low-oxidizer environment some fuel vapor is still left behind 
upon complete droplet gasification, as shown in Fig. 4 . 

4 . Finite-rate kinetics and the flame structure 

It is perhaps somewhat remarkable that because of the basic diffusive nature of droplet 
combustion, much can be learned by making the flame-sheet approximation without having to 
consider the detailed reaction kinetics. There are, however, certain phenomena in which 
the chemical reaction rate can be of the same order as the diffusion rate lit parts of the 
flow field. They can only be described by allowing for finite-rate kinetics and thereby 
resolving the flame structure. Important examples are the determination of the ignition/ 
extinction limits and the combustion of monopropellant droplets. 

To resolve the flame structure, it is necessary to solve Eq. (13) with the reaction 
term, which is governed by the Darokohler number D and the activation energy E. For reac- 
tions of interest to combustion E usually assumes a large value, of the order of 20 to 100 
kcal/mole. Thus the Arrhenius factor exp(-E/R°T) is a very sensitive function of temper- 
ature. It assumes the maximum value at the location of the maximum temperature in the flow 

field. However, by moving slightly away from this location, the slight decrease in the 
temperature can significantly reduce the magnitude of the Arrhenius factor and thereby 
effectively freeze the chemical reaction. 

Thus the flow field can be considered to consist of a narrow reaction region separating 
two broad chemically-frozen regions in which fuel and oxidizer diffuse toward each other. 
Furthermore, since the gradients of the flow properties are very steep within the narrow 
reaction region, diffusion dominates over convection and therefore balances the reaction 
term. On the other hand in the frozen regions diffusion balances convection. Finally, 
since reaction proceeds at a finite rate, the reactants cannot be completely consumed 
within the flame and therefore may leak through it. This leakage reoresents incomplete 

fuel utilization and in severe cases can lead to extinction. The general temperature and 

concentration profiles are shown in Fig. 2b. 

A solution of Eq. (13) usually yields an S-shaped multi-valued curve when a relevant 
parameter, say m, is plotted versus D (Fig. 5). This S-curve can be interpreted as follows. 

For D » 0, the flow is chemically frozen everywhere and we retrieve the pure vaporiza- 
tion state. By gradually increasing D along the lower branch, weak chemical reaction is 
possible leading to slight increases in m. However, with continuous increase in D a value 
Dj will be reached beyond which the only possible solution lies on the upper branch. 

Thus at Di the system changes abruptly from the lower to the upper branch, which repre- 
sents all the possible burning states. In particular as D ♦ » we retrieve the flame-sheet 
solution. Therefore this lower turning point can be identified as the ignition state and 
Di an ignition Damkohler number. Similarly if we continuously decrease D along the upper 
branch, then beyond Dg steady burning ceases to be possible and only the lower branch has 
solutions. Therefore the upper turning point can be identified as the extinction state and 
Dg an extinction Damkohler number. 

Physically, the non-existence of solutions implies that the chemical reaction rate can- 
not balance the heat transport rate. Thus for the lower branch, beyond Dj, the chemical 
heat is generated so fast in the reaction region close to the hot, oxidizing, ambiance that 
they cannot be transported away in a steady manner. This leads to an increase in the reac- 
tion zone temperature and thereby further increase in the reaction rate, culminating in the 
runaway event of ignition. This feed-back mechanism proceeds extremely rapidly because of 
the highly temperature-sensitive Arrhenius kinetics. 

Similarly, for the upper branch, the effects of the finite reaction rate are that the 
flame is broadened and both fuel and oxidizer can now leak through the flame zone. Thus 
with excessive leakage the chemical heat release cannot keep up with the heat transported 
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away from the flame zone, leading to a precipitous drop in the flame temperature and 
thereby extinction. 

Equation (13) has been formally solved using the newly-developed technique of large 
activation energy asymptotic technique. The diffusive-reactive and diffusive-con- 
vective regions are separately solved and matched, using the ratio of the thermal energy in 
the reaction region to the activation energy as the small parameter of expansion. The 
S-curve is reproduced and explicit expressions derived for the ignition and extinction 
states. These criteria are useful in assessing the combustibility of a given droplet- 
oxidizer system. 

Finite-rate kinetics is also essential in the study of monopropellant droplet combus- 
tion. Unlike conventional hydrocarbon fuels which burn through diffusion flames, a mono- 
propellant contains both fuel and oxidizer. Thus after gasification the outwardly-trans- 
ported gas is a premixed combustible, which will eventually burn as a premixed flame. 

While a detailed discussion of the structure of premixed flames is outside the scope of 
this paper, it may be noted that studies* have shown that for large droplets the burning 
resembles that of one-dimensional premixed flame, with 

d(d ) 

- constant , 

dt 

while for small droplets the burning resembles that of the conventional fuel droplet, with 

d(d^) 

• constant . 

dt 


5 . Droplet dynamics 

The heat and mass transfer processes of the spherically-symmetric droplet combustion 
discussed so far are closely influenced by the droplet dynamics, and vice versa. As 
mentioned previously, while increasing non-radial convection generally increases the 
gas-phase transport rates and thereby the burning rate, there exists an upper limit in the 
convection intensity at which the flame can be blown off, that is extinguished. Various 
semi-empirical correlations have been proposed for the burning rate augmentation due to 
convection, such as^ 

K = M + 0.278 Re^^^Pr^^^ , tn(l+B) 

•^s.s. (l + 1.232/RePr4/3]l/2 B 

where Kg_g^ is the spherically-symmetric result. Systematic theoretical and experimental 
investigation on the blow-off of the droplet flame have also been conducted. In particular 
criteria governing the blow-off limits have been determined.*' 

Combustion can influence the droplet drag and therefore its dynamics in tnree essential 
ways. The outward mass flux at the droplet surface reduces friction drag but increases 
pressure drag because of early separation. The net result is as yet unclear. The exis- 
tence of the interfacial velocity, especially after the droplet has been substantially 
heated up, may delay separation and hence reduce both the friction and pressure drag. 
Finally, the significant elevation of the temperature in both the gas and liquid phases, 
and the increase in the gas density because of the presence of the high-molecular-weight 
fuel vapor, can all greatly influence the fluid properties and therefore the drag coeffi- 
cient. The above factors have not been systematically investigated. 

In the section on multicomponent fuels we shall further show that the existence of 
internal motion may qualitatively influence the combustion characteristics of multicom- 
ponent droplets. 


6 . Near- and super-critical combustion 

Because of the enhanced comoustion efficiency and intensity under high pressure, 
internal combustion engines operate under elevated pressures frequently in excess of the 
thermodynamic critical pressure of the liquid fuel. For example, while the critical pres- 
sures of diesel fuels are of the order of 20 atm, th » pressure within diesel engines can 
range from 40 atm at the end of the compression stroke to twice that value at the peak com- 
bustion pressure. While these values seem to indicate that the droplet should undergo 
supercritical combustion early in its lifetime, considerations of the finite droplet 
heating rate and the dissolution of the permanent gases, which can raise the critical 
pressures of the fuel to the range of 40 to 60 atm, render it possible that the droplet 
actually attains criticality after it is mostly vaporized. 


At elevated pressures various aspects o£ low-pressure droplet combustion need to be 
revised. The elevation of the liquid boiling point lengthens the droplet heating period 
and thereby its lifetime. The overall chemical reaction rate generally increases with 
pressure, hence making extinction more difficult but favoring the formation of such pol- 
lutants as NOx and soot. As the critical state is approached, gas-phase compressibility 
becomes important, while the reduced surface tension enhances internal circulation and 
favors droplet deformation and breakup. The gas-phase quasi-steady assumption breaks down 
as the gas and liquid densities become comparable. Finally, at the critical state the dis- 
tinction between gas and liquid vanishes and the phenomena of interest cease to belong to 
the area of droplet combustion. 

7 . Multicomponent droplet combustion 

Much of the earlier studies on droplet combustion used pure fuels. However, recent 
developments in engine design and fuel formulation indicate that multicomponent effects 
will become progressively more important in the utilization of liquid fuels. Combustion 
processes within engines will be more tightly controlled to further improve efficiency and 
reduce emissions. The synthetic fuels derived from coal, tar sand, and oil shale will have 
more complex composition as well as higher and wider boiling point ranges. There also 
exists considerable interest in the utilization of such hybrid fuels as water/oil emul- 
sions, alcohol/oil solutions and emulsions, and coal-oil mixtures. The widely different 
physical and chemical properties of the constituents of these hybrid fuels necessitates 
consideration of multicomponent effects in an essential way. 

To understand heterogeneous multicomponent fuel combustion, either as a droplet or in 
some other forms (e.g. pool burring) , the following three factors have to be considered. 

(i) The relative concentrations and volatilities of the liquid constituents, as would be 
expected. 

(ii) The miscibility of the liquid constituents. This controls the phase change charac- 
teristics. For example the partial vapor pressure of a miscible fuel blend can be 
greatly modified when it is emulsified with evv>n a small quantity of water and 
changes into an immiscible mixture. 

(iii) The intensity of motion within the liquid. This ,'nfluences the rate with which the 
liquid components can be brought to the surface at which gasification takes place. 

The third point requires further amplification. First, it is obvious that no matter how 
volatile a liquid element is, it cannot gasify unless it is exposed at the droplet surface 
either through the passive mode of surface regression, or the active modes of diffusion and 
internal circulation. However, liquid-phase mass diffusion is an extremely slow process; 
its rate being one to two orders slower than that of surface regression. Therefore with it 
being the dominant active mode of transport, it is reasonable to expect that both the vola- 
tile and non-volatile liquid elements i.i the core of the droplet will be "trapped" during 
most of the droplet lifetime. Under this situation the relative volatilities of the indi- 
vidual components obviously cannot be the dominating factors in effecting gasification. On 
the other hand, in the presence of internal circulation, liquid-phase mass transport is 
facilitated such that the relative volatilities exert much stronger influence on the indi- 
vidual gasification rates. 

Thus the combustion characteristics of a multicomponent droplet can be discussed in the 
two extreme rat%p of internal mixing, namely (a) a diffusion limit in which there is no 
motion within the droplet interior such that diffusion is the only active transport mode, 
and (b) a (somewhat artificial) distillation limit in which it is assumed that mixing 
occurs so fast that the states within the droplet interior are perpetually uniformized. 

Let us first discuss the combustion of miscible fuel blends in these two limits. 

Because of diffusional resistance, it is reasonable to expect that the composition of 
the droplet inner core is practically not influenced by the surface gasification and there- 
fore remains close to its initial value. However, since the volatile component has a rela- 
tively lower concentration at the surface, therefore there must exist a thin surface layer 
through which the composition rapidly changes’^* as shown in Fig. 6. Furthermore, 
overall »ass conservation requires that the fractional gasification rate of the ith species 
must be equal to its initial mass fraction. Therefore it is reasonable to expect that 
after the initial droplet heating period, a d^-Law behavior should prevail for the droplet 
size history, as shown in Fig. 7. 

The diffusional stratification of the concentration can lead to the occurrence of an 
interesting disruptive combustion phenomenon termed micro-explosion.^* It is reasoned that 
since the droplet surface has a higher concentration of the non-volatile, high-boillng- 
point, component, and since the droplet temperature is controlled by its composition at the 
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surface, therefore the droplet can be heated to a high temperature. However, the droplet 
inner core is concentrated with the volatile, low-boiling-point, component. It is there- 
fore possible that the inner core mixture can be heated to its limit of superheat and 
homogeneously nucleates. The subsequent Internal pressure build-up is tremendous and 
causes the droplet to explode violently. 

It has also been predicted that micro-explosion can occur only for an optimum range of 
mixture concentration; the non-volatile component is needed to drive up the droplet temper- 
ature while the volatile component is needed to facilitate homogeneous nucleation. ^ ’ 
Furthermore, since the boiling point increases significantly with pressure while the homo- 
geneous nucleation temperature is extremely insensitive to pressure variations, at least 
for pressures not close to critical, it is also suggested*^ that increasing pressure 
enhances the occurrence of micro-explosion. This is an attractive consideration for appli- 
cation in internal combustion engines. 

The occurrence of micro-explosion and its compositional dependence have been experimen- 
tally verified^”, although the pressure dependence has not been experimentally investi- 
gated. 

The possible occurrence of micro-explosion offers interesting potential in the flexibil- 
ity with fuel atomization. That is, fine atomization need not be the primary concern with 
designing the spraying process. In fact, it may be advantageous to have somewhat larger 
droplets which possess sufficient inertial for penetration into the combustor interior in 
order to achieve optimum charge distribution. Upon penetration, rapid gasification can 
then be effected through micro-explosion. 

In the distillation limit the perpetual uniformization of the droplet composition 
implies that the volatility-differentials should be the controlling factor. The gasifi- 
cation of the components should occur approximately sequentially according to their rela- 
tive volatilities. Thus for components with sufficiently different volatilities a plot of 
dg^ versus time should yield approximately linear segments of different slopes separated 
by sharp transition periods, during which the droplet temperature Increases (Fig. 7). Such 
behavior has been experimentally verified. The droplet lifetime may or may not exceed that 
of the diffusion limit. It is also obvious that in the absence of diffusional stratifica- 
tion micro-explosion will not occur in this limit. Thus enhanced mixing, say in the form 
of internal circulation generated by external convection, tends to inhibit micro-explosion. 

It may also be noted that the distillation limit is an artificial one. Recently more 
rigorous analyses* allowing for internal circulation show that the distillation limit can 
never be approached and that diffusional stratif icaton prevails even under intense external 
convection. On the other hand experiments on droplet combustion subjected to either 
natural and/or forced convection do show distillation-like behavior even with mild convec- 
tive intensity.^ As of now no satisfactory explanation exists for this disagreement. 

Next we discuss the combustion of water-in-oil (W/0) emulsions, which have been tested 
in a variety of combustors including the diesel engine, the gas turbine, and furnaces and 
boilers. ’ The results generally indicate a reduction in soot and NO^ emissions while 
those of CO and unburned hydrocarbons are somewhat increased. Change in the combustor 
efficiency on the basis of oil consumption is less clear; slight improvements and degrada- 
tion have been reported depending on the state and type of the combustor, the testing pro- 
cedure, and the properties of the emulsions used. 

The main differences between an W/0 emulsion and a miscible fuel blend is the water 

micro-droplets do not diffuse and that the immiscibility property can significantly influ- 
ence the phase equilibrium relation at the droplet surface. For example under equilibrium 
vaporization, which implicitly requires a distillation mode of internal transport, each 
immiscible component will vaporize independent of the presence and quantity of the other 
component. This implies that the boiling point of the emulsion is limited by the lower of 
the boiling points of these components. Thus the attainable temperature and the associated 
vapor pressure of a high-boiling-point oil can be greatly suppressed by adding only a small 
quantity of water. Nicro-explosion obviously is not possible. On the other hand, in the 
absence of significant internal motion, it has been found^®'^* that micro-explosion again 
occurs and is actually more violent than that for a miscible fuel blend. This is reasona- 
ble because of the large amount of the superheated mass available within the water micro- 
droplets. 

Finally we discuss the combustion of coal-oil-mixtures (COM), which offers a promising 
technology for coal utilization by mixing finely crushed coal in oil and then directly burn 
the coal-oil mixture in such stationary burners as furnaces and boilers. The concept is 
attractive in that it enables the direct substitution of oil by coal, that the crushing and 
mixing processes involved in its preparation are more economical than the energy-expensive 
coal-liquefaction and gasification processes, and that the mixture is still pumpable and 


hence can be readily used in the conventional oil-fired combustors with minimum hardware 
modification. 

The combustion characteristics of COM droplets again depend intimately on the extent of 
internal circulation. Internal bubbling and disruptive combustion should be facilitated 
because, with coal particles serving as potential heterogeneous nucleation sites, the 
droplet does not need to reach the high temperature required to initiate homogeneous 
nucleation. However, the intensity of explosion will be milder because of the reduced 
amount of superheat available. 

Experimentally, sporadic, mild fragmentations have been observed for a suspended 
coal/diesel droplet burning in the quiescent, normal atmosphere. The fragmentation 
becomes explosive with the addition of less than SI water. However, these fragmentations 
can be completely suppressed by increasing the extent of external convection produced 
either by upward blowing of hot gas or simply letting the CON droplet fall tnrough a hot 
furnace. Therefore in the absence of fragmentation, all the coal particles originally 
present in the CON droplet are concentrated into a large coal particle agglomerate^^, which 
requires a long burning time and is expected to be a serious problem with COM utilization 
in combustors with short residence times. 

8 . Droplet interaction 

In the practical situation within a spray, any given droplet is surrounded by and 
therefore interacts with the rest of the droplet ensemble. Since the characteristic 
droplet separation distance can be of the same order as the flame size^, the intensity of 
interaction can be quite strong. In particular, for sufficiently short separation 
distances, the droplets will burn as a merged flame instead of individually. 

To study interaction effects, it is necessary to solve Eqs. (3) and (10), together with 
the momentum conservation equation, for the threo-dimensional flow configuration. If we 
make the thln-flame approximation and also assume weak convection, tlten the analysis 
simplifies to that of the Laplace Equation, or which many solution procedures exist. In 
the presence of convection, special solutions^^"^® have also been obtained. 

Experimentally it has been demonstrated^^ that interaction reduces the droplet burning 
rate because of the competition for onygen. d^-Law behavior obviously does not hold 
because of droplet regression and therefore the continuous variation of the droplet surface 
separation distance. The presence of buoyancy, however, can significantly increase the 
burning rate tecause of the synergistically-enhanced convective transport between 
droplets. 

The case of strongly convective situation is of special interest because of the 
influence of wake regions on neighboring droplets. Another area cf active research on 
droplet interaction is the Leidenfrost Phenomenon^®, in which the droplet gasifies over a 
hot surface and exhibits a non-monotonic behavior in its lifetime with increasing surface 
temperature. 


9 . Concluding remarks 

Presently the fundamental mechanism governing low-pressure, spherically-symmetric, 
single-component droplet combustion can be considered to be reasonably well understood. 
However, quantitative prediction of the burning characteristics is still inadequate, 
primarily as a result of transport property variations. This becomes serious during 
intense burning because the elevated flame temperature not only causes dissociation, but 
also introduces greater temperature and compositional variations in the transport 
properties. It would be cf interest to conduct a numerical simulation of such a process, 
using realistic transport properties but simplified reaction schemes, and compare the 
predictions with the existing experimental data. Caution should be exercised in assessing 
the initial conditions which affect the droplet heating and fuel vapor accumulation 
processes. 

The dynamics of a vaporizing droplet, with various intensities of interfacial mass flux, 
and the coupled problem of droplet vaporization under external convection, have not been 
adequately explored. Systematic studies assuming low, intermediate, and moderately high 
Reynolds number flows are all needed. 

A detailed study of the droplet internal motion and transfer processes could also 
resolve the issue that while theoretical analysis on convective multicomponent droplet 
combustion shows non-uniform droplet concentration persists throughout the droplet 
lifetime, experimental results seem to indicate a batch distillation mode with the 
concentration perpetually uniformized. 
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The phenomenon of micro-exploalon has potential to improve charge preparation and 
deserves further study. For miscible fuel blends the dependence of micro-explosion on the 
fuel composition and internal motion is sensitive and should be characterised. For W/0 
emulsions, emulsion stability, the optimum water droplet sizes, and the effect of pressure 
are the important problems. 

High-pressure near-critical and super-Ckitical combustion is a fertile area of 
research. The primary question here is whether the droplets can reach criticality before 
they have been substantially gasified. Many of the assumptions in modeling low-prersure 
droplet combustion break down near the critical state. The attainment of micro-explosion 
and coal-particle agglomeration also depends on whether criticality is reached before these 
events take place. 

In view of the recent interest to conduct combustion experiments in the gravity-free 
Space-Lab environment, a candidate experiment is that of the combustion and micro-explosion 
of multicomponent droplets. This experiment cannot be easily and perfectly conducted on 
earth using conventional techniques.^ For example the suspension technique cannot be used 
because the suspension fibers can serve as heterogeneous nucleation sites to induce arti- 
ficial micro-explosion. Freely-falling droplets experience forced convection and therefore 
internal circulation, whose intensity also continuously varies because of the changes in 
the droplet velocity and size. Finally, the presence of natural convection also induces 
internal motion. Thus only in a gravity-free environment can a stationary and convection- 
free experiment involving an unsuspended droplet be conducted. 

Recent advances in computational techniques also allow the possibility of studying cer- 
tain aspects of droplet combustion which are not easily amenable to analytical solutions. 
One example is the effects due to variable transport property values as just mentioned. 

The influence is especially strong in the liquid phase. For convective combustion numeri- 
cal solution is also useful because of the lack of spherical symmetry, the presence of the 
complex flows, and the fact that relevant Reynolds numbers frequently are in the analyti- 
cally difficult range of 1 to 100. During high pressure combustion much of the conven- 
tional assumptions describing droplet combustion break down (e.g. gas-phase quasi- 
steadiness and compressibility) such that numerical solutions again may be necessary. 
Finally, studies on droplet interaction may also find numerical approach useful because of 
the lack of symmetry in the problem. 
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spherically-symmetr' droplet combustion. droplet combustion without flow separation. 


GAS-PHASe 



thin-flame 

REGION 


Figure Ic. Schematic showing convective 
droplet combustion with flow separacion. 
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Figure 2a. Temperature and concentration 
profiles during droplet combustion 
in the flame-sheet limit. 
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Figure 2b. Temperature and concentration 
profiles during drop.\et combustion 
with finite rate irinatics. 
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Figure 3. Experimental results on 
the temporal variations of 
diameter-squared, ds^, and 
flame-front stand-off ratio , 
ff=rf/rs» during spherically- 
symmetric combustion. 


Figure 4. Fractional amount of fuel 
reacted as function of fractional 
amount vaporized, illustrating the 
importance of fuel vapor accumulation. 



Figure 5. S-shaped response curve representing all possible solutions 
of quasi-steady droplet burning with finite reaction rate. 
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Figure 6. Droplet concentration profile 
for a bi-component droplet with initially 
equal concentrations. 



Figure 7. Temporal variation of the 
droplet surface area in the two limits 
of internal mass transport. 
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